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An Etude in non-linear Dyson-Schwinger Equations* 
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We show how to use the Hopf algebra structure of quantum field theory to derive nonperturbative results for 
the short-distance singular sector of a renormalizable quantum field theory in a simple but generic example. We 
discuss renormalized Green functions Gn(a, L) in such circumstances which depend on a single scale L — \nq 2 //j, 2 
and start from an expansion in the scale Gh(o, L) = 1 + Ylk 1k{ct)L k . We derive recursion relations between the 
7^ which make full use of the renormalization group. We then show how to determine the Green function by the 
use of a Mellin transform on suitable integral kernels. We exhibit our approach in an example for which we find 
a functional equation relating weak and strong coupling expansions. 
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1. The structure of Green functions 

Our interest is the high energy sector of a renor- 
malizable quantum field theory. We want to ex- 
plore consequences of the underlying Hopf alge- 
bra structure to obtain non-perturbative results, 
using the self-similarity of Green functions. 

Our approach is based on the Dyson-Schwinger 
equations for one-particle irreducible renormal- 
ized Green functions. As we want to obtain non- 
perturbative results the choice of a renormaliza- 
tion condition for us means simply the choice of 
a boundary condition for the full Green function. 

In this short paper we want to exhibit the basic 
idea underlying such a program. We focus on the 
question how to treat the non-linearity of Dyson- 
Schwinger equations systematically. We assume 
we work in a renormalizable quantum field theory 
which provides a finite set 1Z of amplitudes which 
need renormalization. 

For a given superficially divergent amplitude 
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r € TZ we let r r be the sum 
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sym(r) 



(1) 



over all 1PI graphs T contributing to that am- 
plitude, with a a loop-counting small parameter. 
Projection onto suitable form factors cj>(r) allows 
the sum to start with one, so that by the appli- 
cation of the Feynman rules 4>(T r ) is the corre- 
sponding structure function and the Lagrangian 
L is given by 



We can then write PQ 

r = i + B + (r r ,Q({r})),, 



(2) 



(3) 



where Q = Q({T l } r ^n) evaluates to an invari- 
ant charge under the Feynman rules and the 
Hochschild one-cocycle 



B r + {T r , Q) = J2 a k B+ r (T r Q k 



(4) 



fc>i 



is a sum of one-cocycles B^ r and Q is a mono- 
mial in the r r . The uniqueness of Q implies 
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the Slavnov-Taylor identities for the renormal- 
ized couplings pQ. 

The B k / themselves are obtained from the 
skeleton graphs 7 of the theory: 



sym(7) 



B1 



(5) 



where the sum is over all Hopf algebra primi- 
tives 7 contributing to the amplitude r at k loops. 
These maps are defined to be closed Hochschild 
one-cocycles on the sub Hopf algebra generated 
by their concatentations and products |1I2| . 

Effectively, (j3Jl reduces the sum (QJ over all 
graphs to a sum over primitive ones, making use 
of the recursive structure of this fixpoint equation 
which determines the sum of graphs which con- 
tribute to a chosen amplitude. The sums involved 
typically reflect the universal law of |3] and will 
be discussed in detail in upcoming work. 

We set 



r = i + £. 



(6) 



and those Cj are the linear generators of a sub- 
Hopf algebra: 



Theorem 1 There exists maps B_^ r , polynomi- 
als P£j in those linear generators and integers s r 
such that 



• $> fc i?£ r (r-Q*), 



(7) 



AB+ r = B k + r ® I + (id ® B k A A, (8) 



Yl PL ® 4- 



3=0 



(9) 



(10) 



which make the system {c r k } into a sub Hopf al- 
gebra H(A, rn, S, e) of the Feynman graph Hopf 
algebra. 

The polynomials P£ ■ are easily determined and 
we refer the reader to [J for details. 

Feynman rules are then defined in accordance 
with the Hochschild cohomology: 



{Bl(h)) ({<?}) 



d 7 ({k},{q})<P(h)({k})- 



(11) 



Here d 7 ({k}, {q}) is a measure determined by the 
primitive 7 which depends on internal loop mo- 
menta {k} and external momenta {q} such that 
the Hopf algebra primitive determines the kernel 
d 1 and hence the Feynman rules by 



0(BI(I))= / d 7 ({fc},M), 



(12) 



and an appropriate insertion of 4>(X)({k}) in (|11|) 
into the integrand provided by d 7 in accordance 
with the pre-Lie structure of graphs is under- 
stood. Similarly we define renormalized Feynman 
rules by iterating in a subtracted kernel 



[d7({fc},M)-dr(W,M)],(13) 



Bl 



where indicates a suitable renormalization 
point. 

There exists a basis of graphs and external 
structures in the Hopf algebra such that 



(14) 



where L = lng 2 //i 2 is a single scale which deter- 
mines the running of the invariant charge. The 
choice of such a basis disentangles internal subdi- 
vergences into divergent contributions which de- 
pend on a single scale and finite contributions 
which determine the set of primitive elements in 
such a base. In this base, short distance singu- 
larities are captured by Green functions which 
are functions of two dimensionless variables a, L, 
with a remarkable duality between these two vari- 
ables first observed in 

In perturbation theory the Feynman rules now 
allow us to write 

G r R (a, L) = fl (r r ) = l + £ a k Mcl)(L). (15) 

k 

We can expand in a different manner 

G r R (a,L) = l + Y / Y k (a)L k , (16) 

k 

and the renormalization group dictates relations 
between the 7^. We work them out in a moment. 

First, we note that in the case of a linear DSE 
0, we get 



d L <KQ)(L) = 0, 



(17) 
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and hence a scaling solution 

G(a,L) = e -^ {a]L 

solves the linear DSE so that 

7i(a) fc 



lk{a) 



k\ 



(18) 



(19) 



To proceed in general we consider the map 
'W =P lin ® •••®P lin A"- 1 (20) 



where Pn n is the projector into the linear span 
of generators of the Hopf algebra. From |H2j we 
have: 

Theorem 2 The linearized coproduct is obtained 
as 

p nJ rr = p ^r r ® p in r r + Pn n Q ® ad a r r , 

where 

PiinQ = ^2s r P lin T r . (21) 

r 

This allows us to understand the iterative struc- 
ture of the next-to. . . leading log expansion 
We define for n > 1 

' ' ' a, A"" 1 , (22) 



cr„ := — m™ 1 (7i (g) • • 
n! v ^ 



and iTi is the residue defined by 
o-i = d L <j) R (S-kY aug ) {L)\ L=0 . 



(23) 



Actually, a n evaluates to the coefficient of the 
L n term in the evaluation of a Hopf algebra ele- 
ment by the renormalized Feynman rules, by the 
scattering type formula 0]. 

We have 

h £ #Un *i{h) = 0, (24) 

so we can use Theorem 2 and, by the above defi- 
nition of 7t(a), 

1 l(a) = a k (T r (a)). (25) 

Projection onto the linear generators delivers 
the desired formula for the expansion in L: 

7fc(") = \ K(a)7fc-i(a) 
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s J Y 1 (a)ad a <yl_ 1 (a) 



(26) 



With the the above choice of basis we can now 
introduce the Mcllin transform 



FJp)= U(Bl(l))[k 2 



(27) 



(with obvious generalizations to the multivariate 
case as studied below for example) and the DSE 
turns into an equation which determines 7[ as we 
will see in a moment, while the further terms in 
the L expansion are determined from H26J) above. 

The Green function also has the usual expan- 
sion in a which is triangular wrt 



j>k 



(28) 



We can hence proceed to work out the recursion 
relations which express the functions 7^ through 
the functions 7J* for k > 1, and turn the Dyson- 
Schwinger equations into an implicit equation 
which allows to determine the sole unknown func- 
tions 7[(a) from the knowledge of the above 
Mellin transforms. A full discussion is given in 
future work. We now exhibit the approach in an 
example. 

2. A simple example 

For concreteness, we consider massless Yukawa 
theory and consider all self-iterations of the one- 
loop massless fermion propagator, with subtrac- 
tions in the momentum scheme at q 2 = ri 2 . Our 
Green function is an inverse propagator with mo- 
mentum q and a function of two variables a and 
L = lnq 2 //i 2 . We ignore radiative corrections at 
the bosonic line and also at vertices, so the set 
1Z has a single element and the superscript r is 
suppressed henceforth. We rederive the results of 
jSj for this case. 

We write the perturbative series for the Dyson- 
Schwinger equation as 



X(a) = I - aB% 



1 



X(a) 



(29) 



where J (I)) provides the one- loop self- 
energy integral to be iterated. Note that upon 
setting X(a) = I — X(a), this is the equation for 
the self-energy X(a) = — P\i n X(a) studied in [S]. 
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With 

Q = 1/X 2 - P lin (Q) = -2X(a), (30) 

we find the linearized coproduct 

P$X(a) = PunX(a) ® (P lin - 2ad a )X{a). (31) 

This is Proposition 1 of and we also get 

Theorem 3 The next-to next-to. . . leading log 
expansion in L is given through the anomalous 
dimension (a) as 

Jk(a) = rTifaX 1 - 2ad a )j k -i(a). (32) 

This is Proposition 2 of 5 . 

We can now work out with ease the recursions 
which express 7^, k > 1 through the Taylor coef- 
ficients of 71, for example 

72 = \ (7? - 2 7 iad a7 i) (33) 

1 

2 

73 = ^(7i(l-2ad a )7i(l-2ad a )7i) (35) 



[ _ 7ia a2 - 471,171,2a 3 H ] 



(34) 



(3a 3 7i 3 1 +•••)> 



(36) 



and so on. 

Such recursions are obtained for any non-linear 
DSE by iterating Theorem 2. Also, we observe 
that we actually only need the cocommutative 
part in the determination of the coproduct as is 
evident from the very cocommutative definition 
(1221) of o~k, k > 1. The non-cocommutative part 
is always of lower degree in L in the obvious fil- 
tration by L. 

It remains to understand how to compute 
71(a). Instead of an explicit analysis of non-linear 
ODEs as in j^j we proceed here by the Mcllin 
transform, as promised. Here, it reads 



F(P) 



r 
P 



k ■ q 



[fc2]l+p(jfe + g)2 



i>0 



In our simple example we have 
F(P) ' 



p(2-p) 



(37) 



(38) 



Let us introduce a short hand notation: 



7-U = J2lk(a)U k 
fe=i 



(39) 



Then, the Dyson-Schwinger equation becomes 
7 - L = a(l + 7 • a- p )- 1 [e- L " - l]F(p)\ p = , (40) 

where we evaluate the rhs at p = after tak- 
ing derivatives. The functional dependence of the 
non-linear DSE on XQ = X^ 1 reflects itself on 
the rhs. 

The only unknown quantities in this equation 
are the Taylor coefficients 71 j which are implic- 
itly defined through the Taylor coefficients of the 
Mellin transform (|38|l above. 

Taking a derivative of l|40|) wrt L and setting L 
to zero allows us to read them off: 



71 = a(l + 7 VMIp=o 
= ar + a[J2h-d-p} k 



(41) 



k>l 



Lfc=l 

so 71,1 = r universally. 



lj£>=0) 



(42) 



In our example Q = X we furthermore find 



7i,2 
7i,3 



= rf 



rfo+r 2 fi, 



(43) 
(44) 



and so on. 

In this manner one confirms the results of [S] 
which serve here as a mere example for a much 
more general approach. Note that working with 
the toy Mellin transform r — fi = 1 reproduces 
the generating functions which counts the graphs 
contributing at each loop order confirming the 
count of Wick contractions in [5], in terms of 
Catalan numbers 

1,1,2,5,14,42,-... 

For the anomalous dimension itself we indeed also 
confirm the series A000699 of 

A000699 = 1,1,4, 27, 248, 2830, • • • , (45) 

from Ea. l|42|l above, in accordance with [H]. 
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Before we finish this paper by discussing a dou- 
ble Mellin transform we mention one particular 
nice feature of this example which we very much 
hope to work out in general in the future. 

3. A functional equation 

It is our hope that eventually a non- 
perturbative renormalized Green function can be 
related to suitably defined ^-functions. The exis- 
tence of a combinatorial Euler product underlying 
the decomposition of graphs into primitive graphs 
is a first hint [5]. 

Here we report on another such hint based on 
the possibility to reformulate the result of 5 such 
that a functional equation is obtained. Inspection 
of the solution in [H] (which is also suggested by 
the functional equation of the complementary er- 
ror function) shows 



E(o,p) 



VWW) x 

exp(p 2 )erfc(p) 
~ / (exp(p 2 )erfc(p)) 4 



(46) 



where a is now the loop counting parameter and 
p is another variable such that with z = e 2L , 



(47) 



Note that on the lhs of l|4tj[) we have a weak cou- 
pling expansion for a, on the rhs we have a strong 
coupling expansion, hence an expansion in 1/a. 

With T = a/(2n), u = (exp(p 2 )erfc(p))- 4 , 
and Z(T,u) = S(o,p) we get a functional equa- 
tion reminiscent of a functional equation for a 
^-function in two variables for the function field 
case |9lll)| for the non-perturbative renormalized 
Green function 



Z(T,u) 



(48) 



The propagator coupling duality of this Green 
function can now be expressed with u — exp(s + 



t), T = exp(-i), and ((s,t) 
as 

C(M) = -C(t,«), 



exp(^)Z(T, U ) 



(49) 



which we report here for motivation to think fur- 
ther about the connection between quantum field 
theory and ^-functions. 

4. The appearance of transcendentals 

How do we continue in the case where we have 
several insertion places? Nothing changes in the 
above derivation apart from the fact that we now 
have to work with our double Mellin transform 
due to the fact that now both propagators ob- 
tain logarithmic corrections. As the primitive 
self-energy B+(T) now obtains corrections at the 
internal fermionic and bosonic line, we are led 
to consider a coupled system including also the 
bosonic propagator. 

If we are hence to consider the system (based 
now on two elements in 1Z, fermion and boson 
self-energies) 

X(a) = l-aB c , 1 ' 



X(a)Y(a) 
= l-aB c + (X{a)Q(a)), 



(50) 



1 



Y{a) = "nw, 

= I - aB h + (Y(a)Q(a)) , (51) 

with Q~ 1 (a) = X 2 (a)Y(a) and 

fltaQ(o) = -2P lin A» - P nn Y(a), (52) 

describing all possible iterations of massless 
fermion and scalar one-loop graphs, the cor- 
responding functions 7^ (a), 7^ (a) are deter- 
mined in terms of the anomalous dimensions 
"fi(a), y[ (a). The latter arc obtained from the 
system of Dyson-Schwinger equations 



7i 



A" 



a(l+7 X -9-p 1 )- 1 (l+7 y -a- P2 )^ 1 



x (pi + p 2 )F a (pi,/9 2 )U=o, (53) 
7 f = a(l+ 7 X -9-p 1 )- 1 (l+7 Jf x 
x(pi +P2)F b (pi,p 2 )\ Pi = - (54) 

We thus have to consider two Mellin transforms 
in two variables each. They read for the fermion 
self-energy B+(T) 



1 

F c (pi,p 2 ) = — 
T 



-I + P2 



d 4 k 



k ■ q 



eEr=i- 2 C(2fc+i)/ fc (pi,P2) 



(2- pi- p 2 ) (1- pi- p 2 )(pi + pi) 



(55) 
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and for the boson self-energy B b + (I) 
F b ( Pl , P2 ) = ^Jd k [{k + q)2]1+pi[k2]1+p2 

_ 1 - (pi + P2) + PlP2 

(2- Pl - p 2 ) 

e ET=l -2C(2fe+l)/ fc ( Pl ,p 2 ) 

X (1-P1-P2) 2 (P1+P2)' (56) 

where traces have been taken to obtain the rele- 
vant structure functions and where fk{pi, P2) are 
the two- variable symmetric polynomials given by 

t^ = t m Z-i)\ MT 1 -'- (57) 

We observe that the appearance of transcenden- 
tals is utterly generated from the presence of a 
second insertion place. 

Eas. (|53H56(I completely determine the anoma- 
lous dimensions in question and together with 
(|26I52|) the two Green functions. To what extent 
this leads to further functional equations is under 
current investigation. 
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